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We study the superfluidity of single component dipolar Fermi gases in three dimensions using a pairing 
fluctuation theory, within the context of BCS-BEC crossover. The transition temperature T c for the dominant 
p z wave superfluidity exhibits a remarkable re-entrant behavior as a function of the pairing strength induced by 
the dipole-dipole interaction (DDI), which leads to an anisotropic pair dispersion. The anisotropy and the long 
range nature of the DDI cause T c to vanish for a narrow range of intermediate interaction strengths, where a pair 
density wave emerges as the ground state. The superfluid density and thermodynamics below T c , along with the 
density profiles in a harmonic trap, are investigated as well. Implications for experiments are discussed. 


PACS numbers: 03.75.Ss,67.85.Lm,74.20.Rp,74.25.Dw 

I. INTRODUCTION 

Recent experimental realization of quantum degenerate 
Fermi gases of magnetic atoms and the rapid progress 
toward creating degenerate polar molecules 0-0] have 
opened a new frontier for exploring novel phases of quantum 
gases, where dipole-dipole interaction (DDI) plays a central 
role. A lot of attentions have been paid to unconventional 
p-wave superfluids [ ATUII in three dimensions (3D) and topo¬ 
logical superfluids [ 12] in two dimensions (2D). The latter has 
been associated with Majorana fermions and can be used for 
topological quantum computation id. Such exotic super¬ 
fluid phases emerge from the long-range DDI with a strong 
anisotropy, which differs from the widely studied contact po¬ 
tential in dilute atomic gases. Moreover, the relative DDI 
strength can be tuned by changing the fermion number den¬ 
sity n (or Fermi wavevector kp) and, in the case of polar 
molecules 0, by varying an external electric field strength. 

Of particular interest is the intermediate pairing strength 
regime, where complex physics beyond the weak coupling 
BCS theory arises and the superfluid transition temperature 
T c is relatively high, making it more practical to access the su¬ 
perfluid phase experimentally. For a contact potential, the en¬ 
tire BCS-Bose-Einstein condensation (BEC) crossover from 
weak to strong coupling regimes has been studied intensively 
in two-component Fermi gases of 6 Li or 40 K. In contrast, such 
a crossover in dipolar Fermi gases, where richer physics may 
arise, is yet to be explored. Existing theoretical studies in this 
aspect mostly focus on the ground state, based on mean field 
treatments [7, 8,[l5l[lg], which are inadequate in addressing 
moderate and strong coupling regimes at finite temperature. 

In this paper, we address the superfluidity and pairing phe¬ 
nomena of single component dipolar Fermi gases in 3D, with 
an emphasis on the finite temperature and interaction effects. 
Built on previous work lO, [18U that has been applied suc¬ 
cessfully to address various BCS-BEC crossover phenom¬ 
ena in two-component Fermi gases with a contact interaction 
I HI 19], here we construct a similar pairing fluctuation theory 
for the superfluidity of fully polarized one-component dipo¬ 
lar fermions (in the £ direction), in which thermally excited 
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pairs naturally give rise to a pseudogap in the fermion excita¬ 
tion spectrum. We find that (i) the DDI leads dominantly to a 
p^-wave superfluid, and the superfluid T c curve exhibits a re¬ 
entrant behavior as a function of the DDI strength; in the inter¬ 
mediate regime of the BCS-BEC crossover, T c vanishes and 
the ground state becomes a pair density wave (PDW), similar 
to the PDW state studied in underdoped high T c superconduc¬ 
tors im m (ii) In the fermionic regime, the temperature 
dependence of superfluid density and low T thermodynamic 
quantities exhibit power laws, as expected but in stark contrast 
to the contact interaction case [22]. (iii) Within a local den¬ 
sity approximation (LDA), the density profile in an isotropic 
harmonic trap exhibits a similar qualitative behavior to its s- 
wave counterpart, despite the different pairing symmetry and 
the anisotropic pair mass. 

The emergence of the PDW state originates from the long 
range nature of the DDI, which essentially put the system in 
the high density regime. The p^-wave symmetry leads fur¬ 
ther to a non-local effect [23] and hence a diverging coher¬ 
ence length in the nodal xy- plane, which makes it difficult for 
the pairs to move in the £ direction, without heavily colliding 
with each other. At certain intermediate pairing strength, the 
interaction energy between pairs may dominate the kinetic en¬ 
ergy, in favor of forming a Wigner-like crystal in the £ direc¬ 
tion. This PDW state may exhibit behaviors of a Bose metal 
13, 25h . with a Bose “surface” for pair excitations at a finite 
pair momentum q z (with q x = q y — 0). The two dimen¬ 
sionality of the pair dispersion in the remaining xy- plane de¬ 
stroys possible long range superfluid order, leading to a metal¬ 
lic ground state with a density wave of Cooper pairs in the £ 
direction. 

II. THEORETICAL FORMALISM 

We consider an ultracold gas of one-component dipolar 
fermions of mass m in unit volume, with dipole moment 
d = dz, fully polarized in the £ direction. We follow the pair¬ 
ing fluctuation theory as described in Ref. [17], with fermion 
energy = k 2 / (2m) — p measured with respect to the chem¬ 
ical potential p. (we take h = kp = 1, as usual). We shall 
write the pairing interaction into an effective separa¬ 
ble form 12611 . i.e., Fk,k' = ppk<£>£/» where g is the pairing 
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Figure 1. (Color online) Feynman diagrams for the pairing fluctua¬ 
tion self-energy Ti pg and T-matrix t. The thin solid, thick solid and 
dashed lines represent the bare propagator Go , dressed propagator G 
and DDI, respectively. 


process of Taylor expansion, as usual. Following Ref. ini. 
Ti P g can be approximated as T, pg (K) « — Go(— AT)|pk| 2 - 
With the odd parity p_k = —Pk, here we have defined the 
pseudogap A pg as 

A P S = - 2 E^^)~ 2Z_1 Ew> w 

Q q 

where b{pc) is the Bose distribution function. This leads to the 
BCS form of the total self-energy, 

£(K) = X SC (K) + V pg (K) = -A 2 Go(-K)\(fi] i \ 2 , (5) 


strength, pk is the symmetry factor with an odd parity and 
will be determined by the DDI. 

Following previous work 112, 27, [28], the fermion self¬ 
energy comes from particle-particle scattering, which leads to 
both an order parameter (below T c ) and a pseudogap. Non- 
condensed pairs are treated on an equal footing with single 
particle propagators. In contrast to the 5-wave singlet pairing 
case lfl7n . an extra exchange diagram has now been retained in 
the self energy, as shown in Fig. |T] Besides the pairing sym¬ 
metry, this exchange diagram is a major difference between 
singlet and triplet pairing. Therefore, we obtain the fermion 
self-energy from noncondensed pairs 

Z pg (K) = 2%p(K) + £“ chan s e (K) 

= E *(Q) G o(Q - 7f)^k-q/2^k-q/2 

Q 7^0 

- E *(Q) G o(Q - K)tp k—q/2<^3q/2—k , (!) 

Q^O 

where t(Q) = 1 /[g~ l + x(Q)], with x(Q) = 
Ek G{K)G q {Q - K) |^k-q/ 2 | 2 , and G 0 (G) the bare (full) 
fermion Green’s function. Below T c , the condensate self¬ 
energy is 

(K) = -A 2 c Go(-K)\v k \ 2 , (2) 

as in BCS theory, with the superfluid order parameter A sc . 
As in Ref. [ 17], we use a four vector notation, K m (iuj n , k), 
Q = q), Eq = T E*E q > etc -> whh Un (fiz) being 

odd (even) Matsubara frequencies. Here p£ is the hermitian 
conjugate of pk. 

We emphasize that the derivation of this theory is indepen¬ 
dent of the concrete form of the pairing interaction, namely, 
it is not essential whether the interaction is 5-wave, p-wave 
or d-wave, short range or long range, provided that one can 
assume a separable potential in the scattering T-matrix & 
In fact, the original zero temperature BCS-BEC crossover by 
Leggett was done with p-wave pairing [29]. 

Due to the anisotropy of the DDI, the pair dispersion ac¬ 
quires an anisotropy as well, in contrast to the short-range 
contact potential case in a two-component Fermi gas. Namely, 
the finite q pair propagator t pg (Q ) can be expanded as 

tpg (Q) = Dq + P^pair T ^17,q) , (3) 

with an effective pair dispersion D q = + 

g 2 /(2M*) and an effective pair chemical potential p pa ir- 
Here the inverse residue Z and the (anisotropic) effective pair 
mass M]_ = M* = M* and M* can be determined in the 


with a total excitation gap A = y A 2 C + A 2 ^. 

As in Ref. lfl7ll . from the Thouless criteria, t 1 (0, 0 ) = 0, 
we have the gap equation 

and the fermion number equation 


™ = E g w = E 

K k 



+ I 1 f(Ek) 


( 7 ) 


where E k = y/££ + A 2 |pk| 2 is the Bogoliubov quasiparticle 
dispersion and f \x) the Fermi distribution function. 

Now we determine the symmetry factor pk from the DDI, 

V d (v) = d 2 = V(r)Y 2) o(0r,^r), (8) 


where the radial part V(r) = —^16^/5 d 2 /r 3 , and the 
angular part Y2 l o(0 r ,<fi r ) is the spherical harmonic Yi mi ( r), 
with 6 r and <j> r the polar and azimuthal angles of r. So 
the DDI breaks SO(3) symmetry and mixes different partial 
waves. Expanding Vk,k' in terms of partial waves, we have 

Vk,k' = T. W T,rn im ’Crn r {k,k')Y lrni { k)^, (*), with 

= (-l)^l^ 2 w^i\k,k'){lmi\Y2o\l'rn v ) 

and wiji (, k , k') = fj° r 2 drji(kr)V ( r)ji> ( k'r ), where ji(kr) 
is the spherical Bessel function. For a single component Fermi 
gas, only odd l and V are allowed, with V = Z, l d= 2. The r -3 
dependence of the DDI leads to a /c-independent wip(k , k). 
Detailed analyses show that the dominant attractive channel in 
Vk,k is l = 1, mi = 0, i.e., the p z wave, where #oo(^ k) < 0 
is the leading order term, with #oo(^ k) ~ 0.1#oo(^ k) being 
the next leading order term. The leading hybridization terms 
with l = 1, V = 3 are repulsive. Therefore here we concen¬ 
trate on the p z -wave channel. 

To remove the ultraviolet divergence in the momentum in¬ 
tegral of the gap equation, caused by the k independence 
of wip(k,k), we regularize the DDI by multiplying a con¬ 
vergence factor T(r/ro), where ro is the typical radius be¬ 
yond which the DDI becomes dominant [ 13011 . We choose 
F(x) = 1 — e~ x (l + x + x 2 / 2), similar to that used in 
Ref. [15] but here the regularized DDI approaches a finite 
value as r —)> 0, as shown in Fig.[2ta). This leads to a modified 

-wave symmetry factor 



1 - 


ln(l + 4p 2 


4?7 2 


s 2 Ov 


( 9 ) 
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Figure 2. (Color online) (a) Radial part, V{r), and regularized V (r) 
of the DDI, in units of ^/167 t/5 d 2 . (b) k dependence of |^k| 2 cal¬ 
culated from the regularized DDI (black solid line) and |(/?k |nsr (blue 
dashed). For comparison, the radial part of a short-range interaction 
induced p-wave |<pk|p, which scales as k 2 in the low energy limit, is 
plotted as well (red dotted). 


where <pk is real, with g = k/ko = kro , and 9 k the polar 
angle of k. Interestingly, the k dependence of this p>k is quan¬ 
titatively very close to a rescaled 5-wave Lorentzian symmetry 
factor used in Ref. El, 


^k| NS R 


1 + (1.55fc/fc 0 ) 


2 ’ 


( 10 ) 


as shown in Fig. 0b). For comparison, we also plot the k 
dependence of a typical p-wave symmetry factor, \pk\ p = 

-— \oio > induced by a short-range interaction EMI, 

[l + {k/k 0 ) z \ z 

for which the partial wave scattering amplitude f l k ~ Vkk ~ 
|<^k| 2 ~ ai k 21 as fe —> 0 so that for l = 1, a\ is the scattering 
volume. In contrast, the behavior of the p z -wave scattering 
amplitude of the DDI is very similar to the short range s-wave 
case, giving rise to a well-defined scattering length rather 
than scattering volume. Indeed, the strict V (r) gives rise to 
a completely k independent scattering amplitude [8, 34], as is 
the ko —y -f-oo limit of Eq. ©. 

Now with ipk given by Eq. © for the DDI, Eqs. ®, © and 
© form a closed set, which can be solved self-consistently 
for T c as a function of the p-wave pairing strength, g = 
—247rD/(5ra), and for gaps below T c as a function of T, 
where D = md 2 /2 is the dipole length. The unitary limit cor¬ 
responds to the critical coupling strength g c = — 187r/ra£;o, at 
which the scattering length diverges, and a bound state starts 
to form, as determined by the Lippmann-Schwinger equation 
Elim 9c 1 = - Z)k l<y 3 k| 2 /(2e k ), with e k = k 2 /(2 to). 
Thus g/g c = Ak^D/lb. In our numerical calculations we 
take ko/kp = 20, corresponding to a dilute case. 


III. NUMERICAL RESULTS AND DISCUSSIONS 

We first present in Fig. [3] the calculated superfluid transi¬ 
tion temperature T c and corresponding g and pseudogap A pg 
at T c as a function of pairing strength, which are obtained by 
setting A sc = 0. For comparison, the mean-field solution 
T c mf is also shown in Fig. 0a) (red dashed curve). In the 
weak coupling regime, T c follows the mean-field BCS result. 
It starts to decrease after it reaches a maximum around unitar- 
ity g/g c = 1, due to the shrinking Fermi surface. Remarkably, 



Figure 3. (Color online) (a) Superfluid transition temperature T c 
(black solid curve), the mean-field T C MF (red dashed curve) and (b) 
chemical potential g(T c ) as a function of g/g c - Shown in the in¬ 
sets are the pseudogap A pg (T c ) and the inverse pair mass m/M*. 
A PDW state emerges where T c shuts off at intermediate coupling 
strength and the inverse pair mass m/M* becomes negative. The 
right insets share the same horizontal axis as the main panels. The 
lower left inset shows schematic pair dispersion in the PDW regime. 
While the inverse mass remains positive in the xy -plane (blue dashed 
line), it becomes negative in the z direction (black solid curve), with 
a minimum at finite q z in Q*. 


it exhibits a re-entrant behavior. For a range of intermediate 
pairing strength, T c shuts off completely, before it recovers 
at stronger couplings, where the system has entered the BEC 
regime and all fermions are paired, with g < 0. With M* 
approaching 2 m and n pa i r = n/ 2, T c approaches the BEC 
asymptote, 0.1377>, from below. The pseudogap at T c in¬ 
creases monotonically with gjg c . 

In order to understand the re-entrant T c behavior, we plot 
the inverse pair masses in the lower inset of Fig. 0b). It 
reveals that, when T c vanishes at the intermediate pairing 
strength, the effective pair mass in the dipole direction, M*, at 
zero momentum becomes negative, so that the pair dispersion 
D q in the ^-direction becomes roton-like El, with a mini¬ 
mum at a finite q z , as shown schematically in the lower left 
inset of Fig. 0 (solid curve). The pair mass in the xy -plane 
remains positive. This corresponds to a pair density wave 
ground state, with a crystallization wavevector q z in the z- 
direction. Similar PDW states were extensively investigated 
in high T c superconductors in the quasi-2D context II2QL 12111. 

We emphasize that the non-monotonic behavior of T c as 
a function of pairing strength , as found in our T-matrix ap¬ 
proach of the pairing fluctuation theory Mi, can be under¬ 
stood on physical grounds, without invoking specific details 
of the theory. Indeed, this approach has been accepted by in¬ 
creasingly more researchers |33,1380. In the weak coupling 
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regime, T c follows the mean-field behavior. As the pairing 
strength increases towards unitarity, the chemical potential de¬ 
creases, leading to a shrinking Fermi surface and thus a de¬ 
creasing density of state (DOS) N(0) oc ^ p. At the same 
time, a pseudogap develops gradually due to strong pairing 
correlations, which causes a further depletion of the DOS at 
the Fermi level. Both these effects cause a reduction of T c , as 
one can naively expect from the BCS formula for T c . Such 
effects will reach their utmost when the Fermi surface disap¬ 
pears completely at p = 0. Therefore, it is natural to have 
a maximum of T c within the fermionic regime. The actual 
position of the maximum depends largely on the range of the 
pairing interaction, and is close to unitarity in the contact po¬ 
tential limit. On the other hand, as the pseudogap develops, 
fermions form pairs. Upon entering the bosonic regime, es¬ 
sentially all fermions are paired. The BEC temperature of 
these pairs increases with the pairing strength, as the pair den¬ 
sity does. This explains why the combined T c exhibits a min¬ 
imum around p = 0. At this point, the effective pair mass 
M* is significantly heavier than 2m, due to the repulsive in¬ 
teraction between pairs. As the pairing strength increases fur¬ 
ther into the BEC regime, the pair size shrinks, and the inter¬ 
pair scattering length decreases, so that M* decreases gradu¬ 
ally towards 2m. As a consequence, the Bose condensation 
temperature T c of the pairs necessarily increases towards its 
BEC asymptote from below. Within a T -matrix approxima¬ 
tion, these arguments are independent of the specific form of 
the pair susceptibility. 

We note that the emergence of the PDW state has to do with 
the long range nature of the DDI, which essentially put the 
system in the high density regime. At the same time, due to 
the p z symmetry, the coherence length £ ^ up/Ak diverges 
in the nodal xy- plane (i.e., k z = 0) so that the order parame¬ 
ter Ak = A(^k exhibits a non-local effect similar to the case 
of a d x 2 _ y 2 -wave superconductor [23]. (Here vp is the Fermi 
velocity). Such a diverging in-plane coherence length makes 
it difficult for the pairs to move in the z direction , without 
heavily colliding with other pairs. At certain intermediate in¬ 
teraction strength, pairing is strong while the pair size is large, 
so that the repulsive interaction between pairs becomes strong. 
Indeed, a careful look at the effective inverse pair mass reveals 
that before entering the PDW state, the pair mass already be¬ 
comes heavy due to strong pair-pair repulsion. Therefore, the 
kinetic energy of the pairs (in the z direction) becomes much 
smaller than the growing potential energy between pairs, in fa¬ 
vor of forming a Wigner-like crystal structure, which is what 
we call the PDW state. Formation of such a crystal structure 
and minimization of the pair dispersion at a finite momentum 
suppress the superfluid T c down to zero. Such a periodic crys¬ 
tal structure of a PDW state can be most directly probed using 
Bragg scattering, similar to the X-ray diffraction of a crystal 
structure of a solid. 

To further test this picture, we plotted in Fig. [4] the T c be¬ 
havior of the finite range p^-wave superfluid, with a pairing 
symmetry factor given by \pk\ p cos 6>k, as a function of pair¬ 
ing strength for representative values of the range of interac¬ 
tion, as given by ko/kp = 2.5 and 1.0. Here kp/ko serves 
as the effective range of interaction, in units of the interpar- 



Figure 4. (Color online) T c behavior of a finite range p z -wave su¬ 
perfluid as a function of g/g c for ko/kF = 2.5 (black solid) and 1.0 
(red dashed line). The pairing symmetry is given by \pk\ p cos 0^. 

tide distance (1 /kp). For a short range, ko/kp = 2.5, the 
crossover is smooth and continuous, similar to a short range 
5-wave case Hi. except for a reduced BEC asymptote. As 
kp/ko increases, more particles are within the range of inter¬ 
action at the same time so that the effective repulsion between 
pairs becomes strong and the pair mass becomes heavy. For 
a larger range, ko/kp = 1, a reentrant behavior of T c ap¬ 
pears, as in the dipole-dipole interaction case. (And PDW 
states emerge where T c vanishes). In fact, such reentrant be¬ 
havior also occurs for 5-wave pairing with a large range of 
interaction [ 2j|. This supports our conclusion that the reen¬ 
trant behavior of T c for a dipolar Fermi gas results from the 
long range nature of the DDI. We emphasize that the reen¬ 
trant behavior is not unique to the DDI, nor is it to the p-wave 
pairing symmetry. 

In the absence of an underlying lattice potential, the PDW 
state in the dipolar Fermi gases is distinct from a Mott state. 
Instead, it may exhibit behaviors of a Bose metal [ 24,'25|]. The 
presence of the PDW manifests a Bose “surface” for pair ex¬ 
citations Hii , whose energy vanishes at a finite momentum q z 
(with q x = q y = 0). While the pair dispersion remains posi¬ 
tive in the xy plane, the two dimensionality destroys the long 
range superfluid order, leading to a metallic ground state with 
a density wave of Cooper pairs in the z direction. The nature 
of the PDW state deserves further systematic investigations 
& 

It should be mentioned that the chemical potential p 
changes sign within the PDW regime. In the fermionic 
regime, there is a line node at k z = 0 on the Fermi surface 
in the p z - wave superfluid order parameter. Once p becomes 
negative, the node disappears and the excitation spectrum 
becomes fully gapped. This may be regarded as a topological 
transition [[29, Q2]. The anisotropy in the pair mass is a conse¬ 
quence of the DDI. We emphasize that the re-entrant behavior 
of T c is robust against changes of fco and independent of the 
regularization scheme, because fco does not modify the long 
range part of the DDI. It is also present in the next leading 
order, f z - wave channel. 
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Figure 5. (Color online) Transport and thermodynamic behavior, (a) 
ris/n and (b) 7 (T)/ 7 (T C ) as a function of T/T c for g/g c — 0.85 
(BCS), 1.0 (unitary), and 1.5 (BEC), and log-log plot of (c) 1 — n s /n 
and(d) 7 (T)/ 7 (T c ) vs T/T c . 


Note that when /i changes sign, the pairing gap A is rather 
large (of the order Ep). There exists an extended range of low 
T < A, where fi, A and M* remain essentially constant, so 
that the PDW state is rather insensitive to T in this temperature 
range. 

For d -wave pairing as in the cuprates, T c vanishes at a lower 
critical doping concentration, for which the calculated effec¬ 
tive pair mass diverges as well. Below this doping concentra¬ 
tion, the pair dispersion acquires a minimum at a finite mo¬ 
mentum, with a negative mass at q = 0. This suggests that the 
PDW in the cuprates and the PDW in the dipolar Fermi gases 
may share the same origin. 

Next we investigate the transport and thermodynamics be¬ 
havior in the superfluid phase. The superfluid density can be 
derived using a linear response theory. Following Ref. d, 
we obtain 


raA?„ 


yE 

3 


1 - 2 /(£ k ) 
2£ k 


+ f(E k ) 


(v k £ k )Vk| 2 - -(v k f2) • (v k |^> k | 2 ) 


,(ii) 


where f'(x) = df{x)/dx. It can be shown that n s ( 0) = n. 
At 0 < T < T c , both Bogoliubov quasiparticles and pair 
excitations contribute to the thermodynamics. This leads to 
the specific heat C v = ]T k E^dpf(E^) + ]T q D q <9 T fr(^ q ). 

Shown in Fig. [5] are the T dependencies of (a) n s and (b) 
7 = C v /T , for g/g c = 0.85, 1.0, and 1.5, corresponding 
to BCS, unitary and BEC regimes, respectively. These two 
quantities are sensitive to the elementary excitation spectrum. 
Due to the line node on the Fermi surface of the p^-wave su¬ 
perfluid, the low energy density of states N(E) is linear in 
E. Therefore, the low T superfluid density and specific heat 
exhibit power laws in contrast to the exponential behavior of 
an 8-wave superfluid. In the BCS regime, both the low tem¬ 
perature normal-fluid density n n /n = 1 — n s /n and 7 (T) 



Figure 6. (Color online) Comparison of density profiles in an 
isotropic harmonic trap at T/Tf — 0.01, 0.15 and 0.25 and pair¬ 
ing strengths g/g c = 0.85 (BCS), 1.0 (unitary) and 1.5 (BEC). Here 
Rtf is the Thomas-Fermi radius and the density n is in units of kp. 


are linear in T, similar to their counterpart in the nodal d- 
wave cuprate superconductors. On the other hand, in the BEC 
regime, pair excitations dominate, so that n n /n ~ ( T/T c ) 3 / 2 
and 7 ~ (T/Tc) 1 / 2 , similar to the short-range 8-wave case. 
At g = g c , both types of excitations coexist, and thus the T 
dependence exhibits a crossover. The power law behaviors are 
best manifested in log-log plots, as slope changes in Fig.Oc) 
and (d). While the qualitative features shown here may be 
easily anticipated, we emphasize that this is the first system¬ 
atic study of the thermodynamic behavior of a superfluid of a 
dipolar Fermi gas throughout the BCS-BEC crossover. 

Finally, we consider the effect of a 3D isotropic harmonic 
trap of frequency uo with a trapping potential K- ap (r) = 
|mcc 2 r 2 . We assume that Ep is large enough to justify the 
use of FDA lEM l42h . Then /i is replaced by /i(r) = fio ~ 
K-ap(r), where the global chemical potential fi o is determined 
by the total fermion number constraint, N = f n(r)d 3 r, 
with local density n(r). Outside the superfluid core, a non¬ 
vanishing /ipair(^) is included so that the gap and the pseu¬ 
dogap equations are extended as t -1 (0,0) = Zfi pa i r and 
A lg = 2Z -1 ]T q 6(D q - /i pair ), respectively. Shown in Fig. [6] 
is the evolution of the density profile from low to high T, 
throughout the BCS-BEC crossover. Despite the anisotropic 
pairing interaction, the density profile remains isotropic un¬ 
der FDA. It broadens with increasing temperature whereas it 
shrinks with increasing DDI strength, similar to its 8-wave 
counterpart with a contact potential Il42ll . The isotropic den¬ 
sity profile partly reflects the fact that (i) the pairing symmetry 
becomes internal degrees of freedom for the fermion pairs and 
(ii) within the FDA, this isotropy comes from the isotropic 
V[ ra p(r). Possible anisotropy in the density profile may occur 
when direct pair-pair interactions beyond the T-matrix level 
are included, without using the FDA. 

Recent studies using Hartree-Fock approximation, 

suggest that the normal state 3D dipolar Fermi gas is subject to 
collapse and phase separation instabilities in the high density 
and strong DDI regime. For the dilute case considered in the 
present work, the Hartree-Fock contribution to the system en¬ 
ergy, proportional to n 2 , is relatively weak. Our calculations 
show that, within the T-matrix approximation, the compress¬ 
ibility for paired superfluid phase at T < T c remains positive 
definite throughout the BCS-BEC crossover, ensuring a stable 
superfluid state. Effects of direct pair-pair interactions beyond 
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the T-matrix approximation will be investigated in a future 
study. 


IV. CONCLUSIONS 

In summary, our study of single-component dipolar Fermi 
gases reveals a re-entrant behavior of a p^-wave superfluid 
transition T c and a PDW state in a range of intermediate DDI 
strength. Such a PDW state as well as the p^-wave super¬ 


fluid phase may be detected using local density measurements, 
Bragg spectroscopy and momentum resolved rf spectroscopy. 
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